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We compare predictions of the quantum loop expansion to (essentially) infinite orders with (es- 
sentially) exact results in a simple quantum mechanical model. We find that there are exponentially 
small corrections to the loop expansion, which cannot be explained by any obvious "instanton" type 
corrections, ft is not the mathematical occurence of exponential corrections, but their seemingly 
lack of any physical origin, which we find surprising and puzzling. 
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INTRODUCTION 

The Feynman path integral formulationjTj is an in- 
tuitive and powerful method of analyzing quantum sys- 
tems. The lowest order approximation can be understood 
in classical terms, with systematic corrections available 
through a "loop expansion", which is essentially an ex- 
pansion in Planck's constant h. The highlight of such ex- 
pansions is probably the recently completed tenth-order 
QED contribution to the electron^ and muon[3j mag- 
netic moments. The convergence of the resulting series 
is not yet an acute issue for QED, but it is of practical 
interest for theories where the real dimensionless expan- 
sion parameter is much greater, as f.i. QCD. In fact, it 
has been known since the argument of Dyson [1] that a 
power series in the fine structure constant cannot be con- 
vergent, due to instability of QED if a = e 2 / (Airevhc) 
changes sign. An expansion in h is not quite the same, 
but a formal change of sign of h also changes the sign of 
a. Hence, one should not expect more than asymptotic 
series, and hope that they may be given well-defined and 
computable meaning through Borel summation[5]-[7]. 

There are also genuine quantum phenomena, like tun- 
neling processes, which can be understood in quasi- 
classical terms. I.e., as classical processes in imaginary 
time, often referred to as instanton corrections [S]. They 
may lead to non-perturbative contributions which be- 
comes exponentially small as h — > 0. In quantum field 
theory there may also be "renormalon" contributions^ 
which obstructs a Borel summation. However, in sim- 
ple models where the latter phenomena do not occur one 
might think that the loop expansion provides a complete 
description of the computed quantity. At least in princi- 
ple. At least we thought so. 

In ordinary quantum mechanics an expansion in h 
should be equivalent to a WKB expansion (although we 
are not aware of any direct proofs of this). The latter 
seems much simpler to carry out to high orders. The 
WKB expansion can be combined with a quantization 
formula Q first written down by Dunham |10j. which to 



our knowledge has proven to be exact in all cases where 
the result can be computed explicitly to all orders [IT]. 
One might get the impression that Q is always exact 
(Dunham do not claim that). At least we thought so. 
Until we discovered otherwise. 

We have analyzed the perhaps simplest model where 
the WKB result cannot be computed explicitly to all or- 
ders. I.e., the eigenvalue problem 

-r +x^ = e n tp, (1) 

for large eigenvalue numbers N . Here the dimension- 
ed quantity h has been scaled out of the equation. We 
have used 5n = {N + \)~ 2 as the real expansion pa- 
rameter. The WKB quantization formula for systems 
like this was derived to 12 th order by Bender et. aL[TTj 
(counting the standard WKB approximation as th or- 
der). We have recently developed code for very high pre- 
cision solutions of Schrodinger equations in one variable 
(and similar ordinary differential equations) [TH [TJ] . In 
reference fT2l we compared the 12 th order approximation 
of with our very-high-precision numerical computa- 
tions for eigenvalue number n — 50 000. We found agree- 
ment to a relative accuracy of 5 x 10 -67 , which is as 
expected of the WKB approximation for this value of n. 
We interpreted this as a verification of the correctness 
of our numerical code, but it does not constitute a very 
stringent test of the loop expansion itself. 

For a more complete investigation of the latter we have 
extended the WKB approximation to order 1 704. This 
allows us to express the eigenvalues of ([!]) as a series 

E N = const 5 N 2/3 *m^» ( 2 ) 

m>0 

with <5jv = (N + \)~ 2 ■ We have managed to construct a 
double integral representation of the sum in with an 
expression for the integrand which is convergent over the 
whole integration range. For small 6n the result of this 
representation is consistent with an "optimal asymptotic 
approximation" of the sum, i.e. summing the series up to 
(but not including) the smallest absolute value. 
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When comparing results of this approximation with 
numerical computations to about 4000 decimals accu- 
racy, we find an intriguing discrepancy. It vanishes expo- 
nentially fast as (5/v — ► + , in a different manner for even 
and odd eigenvalues. But it is in both cases significantly 
larger than the expected uncertainty in our evaluated 
WKB result. 

The main lesson of our investigation is that the asymp- 
totic series from even very simple model calculations may 
fail to provide results which are "complete" , in the sense 
that they have an accuracy of the same magnitude as the 
accuracy inferred from the optimal asymptotic approxi- 
mation or a Borel summation of the series. 

In the remainder of this letter we present some details 
of our computations and results. 



WKB AND DUNHAM FORMULAS 

For the WKB approximation we formally change equa- 
tion to — e 2 i]f+x i tp — Etp, write ip = e s , and expand 
S = e- 1 J2 n>0 e"S„ to obtain (with Q = V - E) 



— Qi 

n 

S'n—l + ^j^n-j = 0' 
3=0 



(3a) 
(3b) 



which can be solved recursively. The Dunham quantiza- 
tion formula reads 



1 

2i 



S'{z)dz = Nit, 



(4) 



where the integral encircles a branch cut between the 
two classical turning points of |l| at z = ±E 1 l i . All 
odd terms beyond n = 1 in the series for S can be writ- 
ten as the derivative of a function which is single- valued 
around the integration contour [TU1 [TT]. and hence does 
not contribute to the quantization condition Q. The 
even terms may also be simplified by adding derivatives 
of single- valued functions. 



EXPLICIT COMPUTATIONS 

For our case, where V = x 4 , the integral in equation Q 
can be reduced to a sum of integrals of the form 
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r(e) = 

2i J (z 4 _ E) k+l/2 



dz 



(5a) 



= {-^) k \ B(i | - k) E-^-\ for k = 61 - 1, 



r(o) _ 

2i.r (z 4_ jB) fc+i/2 



dz 



(5b) 



M) fc ;|,B(!, i -k)E 1/4 - k , forfc = 6^ + 2. 



Here B(a, b) is the Beta function. These integrals must be 
multiplied by factors {-l) e p^ E 3e and {-\Y p { ° ] E u+1 
respectively, where the pis are p ositive rational numbers 

(e 

2> Pi 

Inserted into M we obtain 



found by solving equation (3b|. F.L, pg = - — 

77 A°) _ 61061 



1768 



Pi 



62 928 ' 



I ff «-l^.-(12f-3)/4 



(6) 



+ 5>l)'B(f , \)qf\^ £"(^+3)/4 = (AT + I) n , 



where and q^ are positive rational numbers. By 
introducing 



[^/B{\,\j\\*E-*l\ 



(7) 



which is a small quantity for large quantum numbers N, 
we can rewrite equation (J6j) as 



1 / oo 



2i 'r 2W e^ \ (8) 



with new coefficients tq = —1, r\ = 1/12tt, and all 
other r m positive numbers growing like m! 2 a™(m + 1)" 
for large m. We have computed these coefficients up to 
m = 852 (selected at random by a license server failure) . 
Empirically they fit the cited behaviour quite well, with 
a r m 0.202 641423 4 and v = — |, and systematic higher 
order correction as a series in m _1 . Obviously the sum 
in equation ^ has zero radius of convergence, but it can 
be turned into a well defined integral by (essentially) a 
twice iterated Borel summation. 

Here we will first invert equation Q to an explicit 
expression for the eigenvalues En. By introducing 5n = 
(N + i)~ 2 we can first express e = En a s a series in Sn, 



(9) 



where the coefficients s m can be computed recursively. 
Their explicit analytic expressions are rational polyno- 
mials in B(j, \ ) and ir^ 1 , which soon become too com- 
plicated for practical use. F.i., 

11 Till II j 5 



S2 



J_ 

' 6ir ' 



S3 



20 736 



144 7T 2 ' 



Instead we have computed s m numerically to about 3 800 
decimals accuracy up to m — 852. Finally we use Q to 
express En by a series in 6n- With the formal expansion 
parameter e = 1, 



E N = [3vr/B(i, i)] 4/3 <5 W 2/3 (l + £ s m 6% 

m>l 

= [37r/B(i,i)] 4/3 S N 2/3 (l+J2t m 6^ 



-2/3 



(10) 
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Borel transformed expansion coefficients B 2 (t,„) (scaled) 



• • • Coefficients for even m 
■ ■ ■ Coefficients for odd m 



(m + i) s / 2 



l<ml 



10 15 20 25 825 830 835 840 845 850 



FIG. 1. Res caled form of the expansion coefficients t m in 
equation ( f 1 for m = 0, . . . , 852. The index v = | is chosen 
as the simplest rational number close to the best fit, after 
which we find at ~ 0.202 64f 423 4 as the best fit to a sequence 
approaching a constant absolute value for large m. 



The first few terms are 



U = 



*2 = - 



648 tt 2 



11 TJ/ 1 1\ 

31 104 7r 4 D V4> 2/ - 



The coefficients s m and i m grow in magnitude in essen- 
tially the same manner as r m , i.e. \t m \ ~ m! 2 a™(?7i + 1)^ 



with at = a r and v = — | (cf. figure [l| . The even and 
odd sequences behave slightly different. Beyond t = 1 
the coefficients t 2 e and £2^+1 have sign (— l) e+1 . 



EXTENDED BOREL SUMMATION 



To make sense of the sum in (|T0| we use the formula 

dyy m e~ a y (11) 



m! 2 = a 2 ( m+1 > 



where a = e 1 ^, with — \n < cf) < \-k. Interchange of sum- 
mation and integration gives a Borel sum (z = xya t S), 



/>oo poo v 

/ dxe~ ax / dye- ay V t r , 

Jo JO ™_n 



t{8) = I dxe 



(12) 



with t m — a 2 ( m+1 ) i m /(m! 2 a™). Our computations in- 
dicate that the sum t(z) = *m zI ™ converge^ for 
|z| < 1. For a = 1 the function i(z) has singularities 



1 We find the equivalent form t(S) = 2 / °° d££ m K (2a^) t(£a t 5) 
to be less convenient for further manipulations. 



where z 2 



-1, with the singular parts behaving like 



(1 + z 2 ) 3 / 2 . In terms of the variable z 2 /(l + z 2 ) the 
points z 2 = — 1 are mapped to 00, and the full integration 
range to the interval [0, 1]. We tried this substitution in 
the hope that the resulting sums for t(z) would converge 
over the full integration range, but discovered additional 
singularites for z 2 ps 4. 

Hence, to avoid integrating through a singularity, one 
must introduce the phase a (or equivalently integrate 
along a different direction in the complex plane) . A con- 
venient choice is a = e 17T ^ 8 , or its complex conjugate. Ac- 
tually, to assure a real result after analytic continuation 
of t(z), one must take the average of these two choices. 



This amounts to taking the real part of the integral ( 12 1 



After our choice of a we separate t{z) into four (in- 
finite) sums, t(z) = J2p=o zP J2e>o^i^+p zii ■ The func- 
tion defined by each infinite sum is singular at z 4 = — 1, 
z 4 ps —16, and probably at infinitely many more points 
on the negative real z -axis. Next we rewrite 



£>0 



- 2J *«H 



= tip) 



t W(z), (13) 



e>o 



and use the computed coefficients t^+p to find equally 
many coefficients t^e+p- A ratio test on the coefficients 



(/-') 



ta+p indicate that the second sum in ( 13 ) con- 
verges for |z 4 /(l + z 4 )| < 1 for all four values of p. This 
provides an expression, 



t{8) = Re 



dx e 



dye 



-ay 



3 

E 

p=0 



z p t^\z), (14) 



where each t^ p \z) is computable by a convergent power 
series in u = z 4 /(l + z 4 ) over the full integration range. 
We have computed i\ p) for t < 212. 



The representation ( 14 ) is not optimal for evaluating 
t(S) for small 5. Instead write e~ ax — —a*-^e~ ax and 
perform a partial integration in x. Repeating this M 



times regenerates the M first terms of the series in (10 1, 



with the remaining coefficients available to construct a 
correction term, 



Af-l 

t(6)= £ t m 5 m + ti^(S). 



(15) 



Here tco2(<5) is an integral similar to t(6) = tcorr(<5). It 
must be computed numerically, but the integral is pro- 
portional to the exactly known prefactor tj^ $ M which 
may be small. A consistency check is that t(S) should 
be independent of M, at least for a range of M -values 
around the minimum of |tcorr(<5)|- As shown in figu re [2 ] 
for the four lowest eigenvalues, the representation (fl5r 
provides results which are independent of M within the 
accuracy of the numerical integration. These results are 
significantly different from the exact eigenvalues. 



(0) 
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Qualitative behaviour of the WKB expansion 
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Logarithmic difference between exact and WKB eigenvalues E N 
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FIG. 2. Behaviour of the WKB expansion of the few lowest 



eigenvalues En- The points show the sum ~^2 r 



tm 



equation ( 15 I. The (green) dashed lines show the full expres- 
sion t(S); it is independent of M to the expected numerical 
accuracy of the correction term tiorl(3)- The red lines show 
the exact eigenvalue, evaluated numerically by our very-high- 
precision routine. Obviously the Borel corrected WKB series 
does not converge towards the exact eigenvalues. 



When investigating a larger range of N we find the rep- 



FIG. 3. Difference between the exact eigenvalues -Ejv, CX act 
(computed numerically to very high precision) and the WKB 
eigenvalues -Bjv.wkb, computed using either the optimal 
asymptotic approximation (OAA) or adding the correction in- 
tegral from Borel summation (Borel) . The results of these two 
methods cannot be distinguished in the figure when N > 1. 
The later is found from equation |7f , with t(<5) computed from 
equation (151 for a range of M-values around 2.2N. The re- 



sult varies little with M, as indicated by the plotted standard 
deviation cj{En)- Hence, the difference between -Ejv, oxact and 
En,wkb is much larger than the uncertainty in -Ejv.wkb due 



resentation (flBl to be consistent with the optimal asymp- to numerical evaluation of the integral for t ( Jf2(S), although 



totic approximation, which is much faster and easier to 
evaluate. For N > 1 the results of the two methods can- 
not be distinguished when compared with the distance to 
the exact eigenvalue. Also, the numerical uncertainty in 



( 15 ) is comparable to the smallest term in ( 10 1 (and even- 
tually worse as N increases). This is shown in figure [3j 
where we plot log l-EW^xact — -Eat.wkbI as function of N, 
together with the expected uncertainties in the evaluated 
values of -E/v,wkb ■ 

For low we find empirically (to exponential accu- 
racy) that 



E 



TV, exact 



E 



TV, WKB 



(-1) 



N 



jN 



(16) 



but, intriguingly, this behaviour is overtaken by larger 
error terms when N > 11 for odd N, and N > 44 for 
even N. We still find £w,exact — ^2JV,wkb to be positive, 
but now the sign of £W+i,exact - E 2 n+i,wkb is 



CONCLUDING REMARKS 

The formula @, with S computed by WKB to all 
orders, does not always provide exact eigenvalues, but 
at best the complete contribution of the WKB approx- 
imation. With hindsight it is clear that there may be 
exponential small corrections: WKB does not predict 



exponentially small as function of N. The correction terms 
look quite simple, with an interesting difference between the 
even and odd eigenvalues. 



backscattering of a forward propagating wave to any or- 
der of approximation when E — V > 0. This is generally 
known to occur in exact calculations. Double backscat- 
tering is likely to contribute an exponentially small cor- 
rection to the left hand side of Q. Also in asymptotic 
analysis a second exponential behavior is said to emerge 
when Stokes lines are crossed, but this statement alone is 
unhelpful for actually computing any exponentially small 
correction. 
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